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Central spin models describe several types of solid state nanostructures which are presently con- 
sidered as possible building blocks of future quantum information processing hardware^. From a 
theoretical point of view, a key issue remains the treatment of the flip-flop terms in the Hamiltonian 
in the presence of a magnetic field. We consider homogeneous hyperfine and exchange coupling con- 
stants (which are different from each other) and systematically study the influence of these terms, 
both as a function of the field strength and the size of the spin baths. We find crucial differences 
between initial states with central spin configurations of high and such of low polarizations. This 
has strong implications with respect to the influence of a magnetic field on the flip-flop terms in 
central spin models of a single and more than one central spin. Furthermore, the dependencies on 
bath size and field differ from those anticipated so far. Our results might open the route for the 
systematic search for more efficient perturbative treatments of central spin problems. 

PACS numbers: 76.20.+q, 76.60.Es, 85.35.Be 



I. INTRODUCTION 

Central spin models are the generic theoretical de- 
scription for several solid state nanostructures which are 
presently under intensive experimental and theoretical 
study in the context of quantum information process- 
ing. Important examples include semiconductor's! and 
carbon nanotube^ quantum dots, phosphorus donors in 
silicon 7 , nitrogen vacancy centers in diamond^!! anc j 
molecular magnets^. The typical Hamiltonian is given 
by 

H = E E r i + E J M + 5 E 5 * (!) 

i—l j — ^ S<.7=1 i—1 

and describes the interaction of N c central spins S% 
with N bath spins Ij characterized by coupling pa- 
rameters Aj with an overall coupling strength A :— 

(1/iVc) J^i^i 12f=i Aj. In semiconductor quantum dots, 
for example, the role of the central spins is played by the 
confined electron spins interacting with the nuclear spins 
of the host material via the hyperfine contact interaction. 
Here the coupling constants A*, are proportional to the 
square modulus of the respective electronic wave function 
at the sites of the nuclear spins and therefore clearly not 
equal to each other ( "inhomogeneous" ) . The parameters 
Jij in the second term of Q account for an exchange 
coupling between the different electron spins, where we 
assume Jjj =: J ex in the following, and the third term de- 
scribes a magnetic field applied to the electron spins. For 
reviews concerning the hyperfine interaction in semicon- 
ductor quantum dots the reader is referred to RefsP^HH]_ 
An important ingredient to the Hamiltonian are 
the so-called flip-flop terms, 

N c N 

i=i j=i 



which arc off-diagonal in the basis with the field direction 
as the quantization axis. Theoretical treatments of ([I]) 
so far have usually distinguished between (i) the case of 
a strong magnetic field, as compared to the overall cou- 
pling strength, B > A, and (ii) the case of a weak mag- 
netic field, B < A. In particular, for the most intensively 
studied situation of a single central spin, N c = 1, the 
flip-flop terms have in case (i) been treated as a pertur- 
bation with A/B being a small parametef^^I], whereas 
in the opposite case (ii) it is commo nly accep ted that 
non-perturbative methods are re quired ! 15 * 23 " ^ However, 
very recently it was shown in RefsP^Ml, again for N c = 1 , 
that surprisingly there is a well-controlled perturbative 
treatment in A(By/N)~ 1 , meaning that, for large enough 
systems, also the case of a weak magnetic field can be 
treated perturbatively. This approach was motivated by 
the statement that the "smallness of the longitudinal spin 
decay is controlled by the parameter A(Q^/N)^ lv (see 
RefpS), where ft denotes the electron spin Zeeman split- 
ting. It is the purpose of the present paper to give a 
systematic and unbiased analysis of the scaling proper- 
ties regarding the flip-flop contributions to the dynamics 
and hence the perturbative regimes. 

II. MODEL AND METHODS 

We want to investigate, in particular, the dependence 
on the bath size N so that we can not make use of 
exact numerical diagonalization. For N c = 1 and if 
Aj = Ay also for arbitrary values of N c , a natu- 
ral alternative woul d be to choose an approach based 
on the Bethe ansatJ^H. This, however, leads to sets 
of algebraic equations which are extremely difficult to 
treat. Therefore, we have to focus on the case of homo- 
geneous couplings, = A/N := A' (see Refs! 24 | 2S | 32 1 ). 
The Hamiltonian Q generally conserves the total spin 
J = S + I, where S := E^i and / := 4 For 

homogeneous couplings it in addition commutes with the 



2 



square of the total bath spin I, 

H, A = \hJ 2 = . 



(3) 



In what follows, we restrict ourselves to the spin length 
Si = Ij = 1/2. We calculate the central spin dynamics 
by decomposing the initial state \a) into eigenstates l^j) 
of the Hamiltonian ([I]) , 



\a) = X^l^)' 



(4) 



and applying the time evolution operator! 24 * 28 ! We will fo- 
cus on initial states with fixed J z quantum number m so 
that only the expectation values of the z-components of 
the spin operators will show non-trivial dynamics. More- 
over, due to the homogeneity of the couplings, the dy- 
namics of the different central spins can be read off from 
each other and we therefore concentrate on the time evo- 
lution (Sf(t)). 

A state which is a simple product of spin states with 
definite z-componcnt is, for vanishing J ex , an eigenstate 
of the Hamiltonian ([TJ) except for the flip-flop terms. 
Thus, for such an initial state all dynamics is due to 
Hff. Therefore, in order to isolate the effect of the flip- 
flop terms we consider initial states of this type, 



|a) = |^...^...fi)®|+... + t .-t), 



NP 



so that 



No 
2 



N 



-N 



D 



Ni 



(5) 



(6) 



Note that for homogeneous couplings the order of the 
spin states within the two subsystems is of no impor- 
tance. 

Since the 2 N dimensional bath Hilbert space is spanned 

by the eigenstates of 1 2 , every product state can be writ- 
ten in terms of these eigenstates: 



ii--4t-t) = EE 



fc=0 {St} 



c^\ N ~k N 



N b D ,{Si}) 

m+Nc/2+NP 



(7) 

Here the quantum numbers {Si} describe a certain 
Clebsch-Gordan decomposition of the bath. Because of 
([3]) , the Hamiltonian |TJ) does not couple states from dif- 
ferent multiplets so that (Sf(t)} decomposes into a sum 
of dynamics on the mul tiplet s given in Q. These contri- 
butions are weig hted bjBBffl 



£ (4 s " 

{s,} 

N°l(N ~ N?)l(N - 



2k + 1) 



(n -k + iy.ki 



(8) 



with k = 0, . . . , -/VjP. In order to compute the dynamics, 
one still needs to perform a diagonalization within the 



2 Nc x 2 Na dimensional Hilbert space of the central spins 
for any value of k in nfh. This diagonalization as well as 
the sum according to 17]) are performed numerically. 

In the following we focus on a weak but finite exchange 
coupling J ex — (1/800)A if not stated otherwise. The 
precise value of J ex is not of significance; similar exchange 
couplings of the same order yield qualitatively the same 
results. However, below we will also briefly comment on 
the special case J ex = 0. The polarization of the central 
spin system or the bath respectively is defined by 



Pc 



Pb 



N r - 2N 1 



N - 2N b D 



N 



(9) 
(10) 



It is well-known that the bath polarization influences the 
cen tral spin dynamics in a way very similar to a magnetic 
fiel d 12 ' 23 ! 28 1 . In the present paper, we will restrict our dis- 
cussions to a very low bath polarization of pb = (1/iV), 
corresponding to N = 2N b D + 1. This is a particularly in- 
teresting special case because any effects of the polariza- 
tion are excluded. However, the results to be presented 
below are clearly not generic with respect to other values 
of Pb- 



III. THE PERTURB ATIVE MEASURES 

In Fig. [l] we give examples of the dynamics for N c — 
1, 2, 3 with the initial states of the central spin system, 
\a c ), given by \a c ) = |M>, IMfT)- In all cases, the 
amplitude of the oscillation is decaying to zero (followed 
by a series of revivals on longer time scales not shown 
in the figures^). The influence of the magnetic field on 
the flip-flop terms manifests itself in two effects so that 
there are two different "perturbative measures" . On the 
one hand, the spin is fixed in its initial direction, in Fig. 
[I] given by (Sf (0)) = —0.5. This means that the magni- 
tude ( "smallness" ) of the spin decay, denoted by n from 
now on, is decreasing with increasing magnetic field. Let 
us further clearify to what extent quantifications of these 
two effects allow to judge about the applicability of per- 
turbative treatments. As explained above, perturbative 
treatments of central spin problems typically consider a 
magnetic field and treat Hff as a small perturbation. 
Clearly, this approach is justified only if the influence of 
Hff on the central spin dynamics is indeed sufficiently 
small. In other words, a perturbative treatment becomes 
the more adequate the stronger Hff is suppressed. 

The quantity \x can be calculated as 



-0.5 -(1/T) / dt(Sm) 



(11) 



which becomes independent of T for T 3> {h/A). The 
measure is illustrated in the first panel of Fig. [TJ On the 
other hand, the decoherence time, denoted by r, as mea- 
sured by the decay of the amplitude of the oscillation, 
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FIG. 1: Spin dynamics for N c = 1,2,3 and N = 401. We 
choose an exchange coupling of J ex — (1/800)A. The mag- 
netic field is fixed to A/B = 4. We consider initial product 
states with \a c ) = |-l|-l|fr) and a very low bath polar- 

ization of pb — (1/N). In all cases, N c — 1, 2, 3, the amplitude 
of the oscillation is decaying to zero (followed by a series of 
revivals on longer time scales not shown). The influence of 
the magnetic field manifests itself in a decrease of the mag- 
nitude of the spin decay and an increase of the decoherence 
time as measured by the decay of the amplitude of the oscil- 
lation. The two quantities, denoted by fi and r, are depicted 
in the first or the second panel respectively. We investigate 
the scaling of the decoherence time by analyzing from which 
time on the amplitude of the oscillation falls under a threshold 
level. The concrete value is of no importance, as long as the 
amplitude is larger than the threshold level before the onset 
of the decay. For high values of p c the quantity r remains 
unaffected by the magnetic field and /j, is the relevant pertur- 
bative measure, whereas for small p c it is r which quantifies 
the influence of the flip-flop terms. 



increases with increasing field strength. In order to cal- 
culate the scaling of r, we fix some adequate "threshold 
level" (Sf) and analyze after which time the amplitude 
falls under this value. The procedure is indicated in the 
second panel of Fig. [I] Note that the concrete value 
of the threshold level is of no importance. It only has 
to be chosen in a way that the amplitude of (Sf(t)} is 
larger than the threshold level before the onset of the 
decay. This procedure has been used already in Ref. 28 . 
Furthermore, a very similar approach has been chosen in 
RefP. 

In Fig. [I] we see that for magnetic fields of identical 
strengths, the value of p, is much smaller for N c — 1 and 
a central spin polarization of p c = 1 than for N c = 2, 3 
with p c = 0, (1/3). Indeed, for large values of p c it is p, 
which adequately describes the influence of the magnetic 
field on the flip-flop terms and for small values it is r. In 
between, both of the measures are of relevance. In the 
present paper we exclude this case and concentrate on 



the three cases shown in Fig. [T] where in case only one 
scale is relevant. 



IV. RESULTS AND DISCUSSION 

Now we come to the central results of the present pa- 
per. In what follows, we investigate the scaling of the 
two measures p, and r with the magnetic field strength 
for a fixed particle number and with the particle number 
for a fixed magnetic field. With respect to the scaling of 
p, we consider N c = 1, 2 and for r we focus on N c = 2, 3. 

As already mentioned above, due to (|3| and Q, the 
dynamics (S[ (t)) decomposes into a sum of dynamics on 
different multiplets. In a first step, it is instructive to 
focus on N c = 1 and to consider only a single term of 
the sum. As to be demonstrated below, this leads to the 
scaling of the measure /i with the magnetic field on a 
fully analytical level. 



A. Magnetic field scaling of /j, for N c = 1 

As we are dealing with only a single central spin in this 
subsection, in what follows we drop the index in (Si(t)). 
Let / denote the quantum number of some multiplet in 
the sum (S z (t)). This corresponds to the Hamiltonian 



H, = A'S ■ I + BS Z . 



(12) 



For fixed m this corresponds to a 2 Na x 2 = 2x2 matrix, 
which can be diagonalized easily. Here it is convenient 
to respresent Hi with respect to the eigenbasis of the 
first term, S ■ I, resulting from the well-known formula 
for coupling a spin of arbitrary length to a spin of length 
5 = 1/2 (see e.g. RefP) 



\I±\,m) 



= c±(m)\tf)\I,m- 
± c T (m)\^)\I,m- 



where 



c±(m) = 
This yields the matrix 



I±m + l/2 
21+ 1 ' 



(13) 



(14) 



/ #f _|_ Bra 

2 *~ x 



Hi 



B 2 _ B 2 m 2 \ 
4 x 2 



IE 

4 



B 2 m 2 



A'(I+1) 



Bra 



(15) 



7 



where we introduced the shorthand notation x = 2J+ 1. 
We denote the components of the Hj eigenstates with re- 
spect to the basis {If) | J, m-1/2), |JJ.) \I, m+l/2)} by ipV ) 



4 



and the corresponding eigenvalues by Diagonalizing 
( 15 1, we get 



with 



V4 
4 2) 



a± 



= (a + c~ (771) + fe + c + (m)) 

= (&+c~(m) — a + c + (m)) 

= (a_c _ (m) + 6_c + (m)j 

= (fo_c~(77i) — a_c + (rn)) 



Here we defined 
and 



= T 



1 

" 4B 2 x 2 z 2 + z 2 _ 

1 + {A'x 2 +iBrn±A'xy) 2 

1 Bz + z_ 



z± 



AB 2 8Bm 



A' 2 



x ± 2m 

x 



(16a) 
(16b) 
(16c) 
(16d) 

(17a) 
(17b) 

(18) 
(19) 



Considering the initial state \a) = |JJ.) \I, m+1/2), it then 
follows for the central spin dynamics: 



(S*(t)) = 



(20) 



(1)|2 



(2) |2 



(2) 1 2 



+V4 V: 



cos 



(^ - E 2 )t 



We denote the measure \i corresponding to ( 12 ) by 



Obviously, the quantities , introduced in (|20|), are 
related to fi^- 1 ' by 



1 



(21) 



Inserting jl7| in the expression for fJQ given in (20 1 and 
performing some extensive algebra (see appendix) , we get 

(J) _ 1 (B/2A') + (m/2) (x 2 - 4m 2 )/4 
^1,2- 4 T y v 2 { ' 



and hence 



m _ Qr 2 -4m 2 )/2 



(23) 



Let us consider i? to be given in units of A (for sim- 
plicity we denote £? = _BA). Then we have y 2 = 
x 2 + 4BN(BN + 2m). As mentioned above, in our anal- 
ysis we focus on initial states with nearly unpolarized 
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FIG. 2: B field and iV scalings of fj, for AT C = 1,2 and 
N = 401,801 or A/B = 8,4, respectively. The exchange 
coupling in the case of N c — 2 is fixed to J ex — (1/800)A 
We consider initial product states with \a c ) = j i}-), i 
corresponding to p c — 1, and a low bath polarization of 
Pb — (1/N). The results are plotted on a double logarith- 
mic scale. We find power laws ~ B~ v with v « 2 and 
~ AT - " with ;/ « 1. In the first case, the exact values are 
given by v = 2.02485,2.09705 (N c = 1,N = 401,801) and 
v = 1.95565, 1.96427 (N c = 2, N = 401, 801). For the N scal- 
ing we have v = 0.989959, 1.01776 (N c = 1, A/B = 8,4) and 
v = 0.925328, 0.97564 (N c = 2, A/B = 8, 4). With respect to 
the magnetic field scaling, the fully analytical result, given in 
(241, is reproduced. 



baths. The measure fi is significant only for highly po- 
larized central spin systems. Hence, in the most impor- 
tant situation of N c <C N it follows BN ^> 2m so that 
y 2 scales like B 2 . Consequently, we have fi^ 1 ' ~ B~ 2 . 
This result is independent of the value of / and, hence, 
we arrive at 



1 

B 2 ' 



(24) 



B. General scaling properties 

Now we evaluate the full dynamics in an almost ana- 
lytical fashion and derive the scalings of fi (for N c = 1,2) 
and t (for N c = 2, 3). As to be demonstrated below, the 
number of central spins has no (direct) influence on the 
result. In Figs. [2] and [3] we plot n and r against the 
magnetic field and the number of bath spins on a dou- 
ble logarithmic scale. We consider two different values 
of N or B respectively for each number of central spins. 
Obviously, this changes the values of the measures, but 
not their scaling properties. For the B field scaling we 
find a simple power law ~ B~ 2 in all cases, which for 
N c = 1 reproduces the fully analytical result presented 
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FIG. 3: B field and TV scalings of r for TV C = 2, 3 and TV = 
401,801 or A/B = 8,4, respectively. The exchange coupling 
is fixed to J ex — (1/800) A We consider initial product states 
with |q c ) = ID-'fl - ), j JJ-JJ-Tr) , corresponding top c = 0, (1/3), and a 
low bath polarization of pt — (1/TV). The results are plotted 
on a double logarithmic scale. We find power laws ~ B v with 
v m 2 and TV" with v ~ 2. In the first case, the exact values 
are given by v = 2.07975,1.99543 (TV C = 2, TV = 401,801) 
and v = 1.97742,1.99555 (TV C = 3, TV = 401,801). For the 
TV scaling we have v = 1.87736,1.97722 (TV C = 2, A/B = 
8,4) and v = 1.77072,1.92503 (TV C = 3, A/B = 8,4). Note 
that an increase of r corresponds to a decrease of the flip-flop 
contributions to the dynamics. 



above. Indeed, already this is much stronger than the 
B _1 scaling ant icipated by the perturbative approaches 
presented so f ar pIH211 2 9 | 30 | j\f te that an increase of r in- 
dicates a decrease of the influence of the flip-flop terms 
on the dynamics. The scaling with the number of baths 
spins turns out to be even more surprising. For [i we find 
r-j TV -1 , whereas for r the influence of Hff scales down 
with ~ TV~ 2 . As mentioned above, for TV C = 2,3 we al- 
ways considered a weak but non-zero exchange coupling 
Jex — (1/800) A The results are generic for J ex ^ 0. 
However, for a zero exchange coupling, J ex = 0, the r 
scaling yields a slightly different result. Here the expo- 
nent in the magnetic field scaling ~ B~ v decreases from 
v = 2 to v = (3/2). We have no explanation for the 
concrete value of 3/2. However, it is not suprising that 
v does not become 2 as in the case TV C = 1, where we 
naturally have J ex = 0. This would be the natural guess 
if for J ex = the Hamiltonian would decompose in a 
sum of independent TV C = 1 models. However, this is not 
the case as the TV C central spins interact with a common 
bath. Here the dynamics of the central spins result from 
the interaction with the spin bath and with each other 
through the spin bath. The c ase of separate baths has 
been investigated in RefsF^Il. 

Consequently, in any case the somewhat surprising ap- 
proach to use A{B\/N)~ 1 as a the small parameter for a 



perturbative treatment, presented in Refs 1221251 turns out 
to even underestimate the flip-flop suppressing character 
of the particle number. Hence, with respect to the par- 
ticularly interesting low field case the perturbative limit 
is not yet achieved. 



V. CONCLUSION 

In conclusion, we have studied the scaling of the influ- 
ence of the flip-flop terms with the magnetic field and the 
number of bath spins in central spin models. In order to 
be able to treat comparatively large systems, we consid- 
ered homogeneous couplings. The flip-flop contribution 
to the dynamics has been isolated by choosing simple 
product initial states. The effect of an applied magnetic 
field manifests itself in the magnitude of the spin decay 
and the decoherence time. For highly polarized central 
spin systems it is the magnitude of the spin decay which 
is the relevant scale, whereas for a low central spin polar- 
ization it remains, to a large extent, unaffected and the 
decoherence time describes the influence of the magnetic 
field. We investigated the scaling of fj, and r for TV C = 1, 2 
and TV C = 2, 3 in different parameter regimes. 

Surprisingly, we found that /i decreases quadratically 
with the magnetic field and linearly with the number of 
bath spins. For TV C = 1 we presented a fully analyti- 
cal derivation. For J ex > the decoherence time shows 
identical scaling properties with respect to the magnetic 
field, whereas if J ex = 0, the behavior slightly changes to 
~ S" 3 / 2 . As a very interesting and unexpected result, it 
turns out that r increases quadratically with the number 
of bath spins, corresponding to a quadratic decrease of 
the flip-flop contributions to the dynamics. Summariz- 
ing, for p c w 1 our results suggest 

A{By/N)- 2 

as the small parameter of a perturbation theory. This 
esse ntially goes along with the approach considered in 
RefsPSGEI. However, for small values of p c we have 

A(BN)- 2 . 

This means that the perturbative treatments of central 
spin models presented so far strongly underestimate the 
influence of both, the magnetic field as well as the num- 
ber of bath spins. It is therefore desirable to search for 
new approaches using the full suppression of the flip-flop 
terms. Although all scaling properties are independent 
TV C , central spin models with more than one central spins 
are particularly interesting with respect to such investiga- 
tions, as here low polarizations of the central spin system 
can be achieved. This leads to a very strong decrease of 
the influence of the flip-flop terms with the number of 
bath spins and hence the possibility to treat extremely 
small magnetic fields. 
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(1 1 2 h in the expression for ^\ given in (|20|) we get 



Appendix 



In what follows we present details on the derivation of 



(22). Inserting the eigensystem (16) of the Hamiltonian 



a {1) 
An, 2 



which can be simplified to 



_(B + A'm) (A'x 2 + ABm ± A'xy) (Ax 2 ± A'xy + 2B (2m + xz 2 _)Y 
4x(A' 2 x 3 ± A'x (±8Bm + A'xy) + 4B (B + B(x - 1) ± Amy)) 2 



(25) 



(/) -(B + A'm) (x + 2m) (A'x 2 + 2Bx ± A'xy)' 

Ml,2 = 



4a; 2 A' 2 y 2 (4Bm + A'x 2 ± A'xy) 



(26) 



On the other hand, the expression (22) can be rewritten 
as 



(7) _ (B + A'm) (B + A'm) 2 
Ml ' 2 ~ T 2Ay + A' 2 y 2 



(27) 



If we now equalize (26) with (27) and multiply by the 



denominators, we get for one side 



- 165 2 m =f SA'Bmy - WA'Bm 2 
=F 2A' 2 x 2 y - 4A'Bx 2 - 4A' 2 x 2 m 

- 2A' 2 xy 2 =F 4A'Bxy =p 4A' 2 mxy 



and for the other 



- A' 2 x 3 T A' 2 x 2 y - 2A'Bx 2 
=F A'Vy - A' 2 xy 2 - 2A'Bxy 

- 2A'Bx 2 =F 2A'Bxy - 4S 2 x 

- 2A' 2 x 2 m =F 2A' 2 xym - AA'Bmx 
=F 2A' 2 xym - 2A' 2 y 2 m =p 4A' Bmy 

- AA'Bmx =F 4A'Bmy — 8B 2 m 



(28) 



(29) 



Inserting ( 18 ) in the terms proportional to y 2 immedi- 
ately shows that (28) and (29) are identical, which yields 



(22) 
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